Magneto-mechanical interplay in spin-polarized point contacts 
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We investigate the interplay between magnetic and structural dynamics in ferromagnetic atomic 
point contacts. In particular, we look at the effect of the atomic relaxation on the energy barrier for 
magnetic domain wall migration and, reversely, at the effect of the magnetic state on the mechan- 
ical forces and structural relaxation. We observe changes of the barrier height due to the atomic 
relaxation up to 200%, suggesting a very strong coupling between the structural and the magnetic 
degrees of freedom. The reverse interplay is weak, i.e. the magnetic state has little effect on the 
structural relaxation at equilibrium or under non-equilibrium, current-carrying conditions. 
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PACS numbers: 75.75.+a, 73.63.Rt, 75.60.Jk, 72.70.+m 

Existing experimental techniques are capable of con- 
structing low-dimensional magnetic constrictions consist- 
ing of a tiny number of atoms and measuring their trans- 
port properties 0, 0, Q I n this rapidly evolving field 
experiments and theoretical simulations arc closely en- 
tangled and the latter can provide important informa- 
tion on the fundamental dynamics at the atomic scale. 
The modeling of atomic-sized ferromagnetic devices un- 
der bias requires the combined description of electron 
transport and of the local magnetization and structural 
dynamics at the atomic level In fact, since the spin- 
polarized current transfers both spin and charge, it exerts 
a twofold effect on the current-carrying structure, inter- 
linking its structural and magnetic degrees of freedom. 

In order to investigate the mutual interplay between 
magnetic and structural properties, we have developed 
a computational scheme for evaluating spin-polarized 
currents and the associated current-induced forces and 
torques. On one hand, we are able to examine the ener- 
getics of magnetization dynamical processes as a function 
of the atomic rearrangements in magnetic point contacts 
under bias. On the other hand, we can ravel out the 
effect of the magnetic configuration itself on the struc- 
tural relaxation. Our main finding is that this interplay 
is strong only in one direction. While the atomic re- 
arrangements can modify drastically the point contact 
magneto-dynamics, the magnetic configuration has little 
effect on the atomic configuration. 

Our computational approach combines the non- 
equilibrium Green's function (NEGF) method for eval- 
uating the charge density and the current 0, Q with a 
description of the magnetization dynamics in terms of 
quasistatic thermally-activated transitions between sta- 
tionary configurations Q. This scheme is general and 
is conceptually transferable to first-principles Hamil- 
tonians, for instance [|J, within density functional the- 
ory. However, it is currently implemented in an empirical 
single-orbital tight-binding (TB) model || with parame- 
ters chosen to simulate the mechanical properties of no- 



ble metals and the basic electronic structure of a generic 
magnetic metal. This has the benefit of being reasonably 
realistic while keeping the computational overheads to a 
minimum. 

The magnetic and structural degrees of freedom both 
enter our model as classical variables. We associate a 
classical magnetic dipole moment (MM) and a set of 
Cartesian coordinates to each atom in the system. The 
magnetic state (MS) is defined by a set of angular coordi- 
nates $ = {4>i}- These represent the angles of the MMs 
of the "live" atoms (those contributing to the dynamics) 
with respect to a given direction. The structural state 
(SS) is defined by a set of spatial coordinates 1Z = {Ri} 
for these atoms. The interplay between MS and SS is 
then investigated by either fixing $ and evolving 1Z or 
by fixing 1Z and evolving <!>. In the first case we study 
how the MS affects the structural relaxation, and in the 
second how the SS modifies the magnetic dynamics. 

Our mixed quantum-classical Hamiltonian for this sys- 
tem reads 



H{*,K) = X)[ fl tf B W + W] c U- + 
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where cj CT and Ci CT are creation and annihilation operators 
for electrons with spin a at atomic site i. H(^,TZ) is 
separated into several quantum and classical terms. The 
first term reads 



(2) 



where £q is the on-site energy for an isolated atom, V? B 
represents an empirical approximation to the second- 
order Coulomb energy at the i-th site 
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where we have defined fik(lZ), n e i = e 2 /47r£ = 14.4 
cVA, Aqj = J2a Pjj ~ (Po)jj is the deviation of the elec- 
tron occupancy of site j (pjj being a diagonal term of 
the density matrix for spin a) from its equilibrium non- 
interacting value pq. The last term in Eq. @ is the 
intersite hopping integral, depending as an inverse power 
q = 4 on the interatomic distance Rij = \Ri — Rj\. Here 
we use e = 7.8680 meV, c = 139.07, a f = 4.08 A, which 
are parameters fitted to elastic properties of bulk gold 
H . The value of the on-site Coulomb repulsion is U = 7 
cV and the band filling is po = 0.7 electrons per atom. 
This non-integer band filling effectively accounts for the 
s-p-d hybridization. 

The structure of our atomic point-contact system is 
presented in Figure ^ It consists of a linear chain of 
three atoms attached to two semi-infinite leads with a 
simple cubic structure and 3x3 atom cross-section. The 
lead magnetizations are polarized along the z-axis and 
they are anti-parallel (AP) to each other. The MMs 
of the atoms in the chain are allowed to rotate in the 
(x — z)-plane, and are described by the set of polar co- 
ordinates <f> = {<fii} (i = 1,2,3) about the z-axis. In 
this setup an abrupt domain wall (DW) is formed in the 
constriction. Note that our model does not include spa- 
tial spin anisotropy, therefore the direction of the spin- 
quantization axis z is arbitrary. In this sense a Neel and 
a Bloch wall arc physically identical within our model. 




FIG. 1: (Color online) Schematic atomic configuration of the 
magnetic point contact investigated in this work. Some dis- 
tances are magnified for clarity. 

The spin-dependent term in the electronic part of Eq. 
is approximated by a classical Heisenberg-type inter- 
action as follows 

V? (*) = -Ja ■ SAj = -o z J - cos (^Sij , (4) 

where a = (0, 0, \cr z ) is the spin of the current-carrying 
electrons, collinear and polarized along the z axis (a z = 
±1) and J > is the exchange integral. We assume 
that the local MMs, associated with the local spins, arc 
constant in magnitude 0, and impose |5,-|=1. 



The remaining terms in Eq. involve only classical 
variables. We have adopted a repulsive pai r potential 
which decays as an inverse power law 0] 

si(n) = \ Y, % m = \ E (|f:) . (5) 

and a Hciscnbcrg-type nearest-neighbors spin-spin inter- 
action between the localized spins 

W{*) = - J -f £ Si ■ Si = - J -f £ cos {4> t - ^ , 

(6) 

where Jdd > is a parameter describing the strength of 
the intersite exchange interaction. Typical values for the 
exchange parameters are J = 1 eV and Jdd = 50 meV, 
which are of the same order of magnitude as those for 
bulk ferromagnetic metals within the s — d model |l3j . 
and those derived from the Curie temperature within the 
Heisenberg model Clearly, the low coordination al- 
ters the values for the exchange parameters 0, though 
this effect was not explicitly included here. However, we 
have repeated the calculation for J from 0.8 eV upto 2 eV 
as well as for other force parameterizations (for Cu ||) 
and found no qualitative differences in the main physical 
results. 

Forces and torques, associated with the classical 
degrees of freedom, are derived from the Hcllmann- 
Feynman theorem, generalized for non-equilibrium sys- 
tems. 0, llOj In our model the corresponding forces 
read[ll| 

F i = -Y, [ 2 (V<-H£ B )Refo,] + A-/,A,/,V,/, ; + Vifitf] , 

(7) 

where the live atoms (i = 0, . . . , 4) are those in the chain 
and their leftmost and rightmost neighbors in the leads 
(i = 0,4). Index j runs over all the atoms in the self- 
consistent region, which apart from the chain includes 
two atomic planes from each lead (see Fig. The 
torques T read 

Ti = -^Sjsin^j-^ [sin(0 4 - + sin(0 i+ i - 4> t )] , 

(8) 

where Si = p\ — p\ { is the on-site spin-polarization (i = 
1,2,3), 0o = and <pi = it are considered frozen and 
aligned along the same direction as in the corresponding 
lead. 

The quantum transport problem is solved by using the 
NEGF method 0, 0, • ^ n analytical expression is used 
for the surface Green's function of the leads and a positive 
external bias V is introduced as a rigid shift of V/2 of 
the on-site energies of the left-hand side electron reservoir 
and of — V/2 for the right-hand side with respect to their 
equilibrium value. The non-integer band-filling of the 
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spin-split s band in the leads generates an asymmetry 
in the transport properties of the spin-up and spin-down 
conduction channels upon left-to-right inversion of the 
system in its AP state. The relative spin-polarization 
P a = (rrf — n^)/(n^ + n^-) of the density-of-states at the 
Fermi level, n CT (£*/), of the left lead in equilibrium is 
about -19% (-17%) for J = 1 eV (2 eV). The net current 
is calculated as in Ref. Q ■ The self-consistent density- 
matrix is used into Eqs. (0 and ijHJl, and both forces and 
torques are determined. 




FIG. 2: (Color online) Displacements (in picometers) of the 
atoms in the chain: (a) from the uniform geometry at V = 0; 
(b) as function of the DW migration reaction-coordinate 02 
at V = 1 V; (c,d) as function of the bias voltage V. See text 
for details. Here J = 1 eV, Jdd = 50 meV. 

The "live" atoms in the constriction are relaxed at a 
given finite bias for a fixed MS. For symmetry reasons the 
atomic relaxation always results in displacements only 
along the longitudinal direction (y-axis). Thus, we de- 
note the relaxed atomic positions at a bias voltage of V 
and a magnetic state during relaxation by yy(0). The 
initial geometry, denoted by y U nU is that of equidistant 
atoms with a nearest- neighbor distance of a = 2.5 A (see 
Fig. [IJ. This is near the equilibrium bond length of a 
periodic ID chain within our model |l2^ . However, such 
a bond length produces a compressive stress in the bulk 
leads, as a result of which the leftmost and rightmost 
atoms are pushed slightly out of the leads and the whole 
chain at equilibrium shrinks by about 2% [see Fig. Eta)]. 

First we investigate how the MS of the constriction af- 
fects its structural relaxation. We relax the live atoms 
in {0^} = (0,0, 7r) or {0^} = (0,7r,7r) configuration, 
which represent two possible spatial positions of the DW 
inside the constriction, as well as for intermediate MSs 
with 2 G [0, 7r] and 0i i3 such that T lj3 = 0. The DW- 



motion induced displacements A?/^ = yv{4>2) — 2/v(007r) 
are monotonic functions of 02 [see Fig. EJb), where 
V = IV"]. They represent nearly rigid translations of 
the whole atomic chain in the direction of the electron 
flow. However, these displacements are very small and 
constitute about 3% of the displacements from the uni- 
form structure Ay uni = yo({0 (1 ' 2) }) - 2/uni [Fig. Ha)]. 
The former depend weakly on the bias [Fig. Eld)]. 

The effect of the bias voltage on the atomic relaxation 
Ayy = 2/v(0) — yo(0) for = {0^ 1,2 '} is demonstrated in 
Figure El^c). The relative current-induced displacements 
are smooth functions of V and promote a tendency to- 
wards dimerization. \Ayy\/a < 0.7% (upto V — 2 V) 
are of the same order of in magnitude to the displace- 
ments Ay^. Interestingly, the former are almost insen- 
sitive to the MS {0' 1 ' 2 -'}, which is a result of the weak 
bias-dependence of Ai/r^n [see Fig. Efd)]. Other MS, 
having multiple abrupt DWs [e.g. (ir, 0, ir)], are found to 
produce an effect of similar magnitude on the atomic re- 
laxation (not shown in this work). This establishes that 
the structural properties of a magnetic nano-device un- 
der bias are to a large extent independent of the magnetic 
state. 

Next we address the converse, that is whether the 
structural relaxation affects the energy barrier for DW 
migration. This is calculated from the torque T2 as the 
MM of atom "2" is quasistatically rotated. The other 
two live MMs are continually relaxed so that their asso- 
ciated torques are kept at zero [?J ■ The net work for this 
rotation is then 

VF(0 2 ) = - r 2 T 2 (0 2 )d0 2 (9) 
Jo 

Our calculated DW migration barriers (DWMB) for a 
few values of J are shown in Figure [3J A key feature of 
these profiles is the tilt of the barrier at any finite bias 
even for a uniform atomic arrangement. This tilt results 
from the fact that in the absence of the special symme- 
try in the density of states that characterize the present 
simple cubic structure in the case of a half-filled band , 
the leading contributions to current-induced forces and 
torques are linear in the bias. This asymmetry results in 
a preferential spatial localization of the DW at a given 
bias: under the present (positive) bias the (0, 0, ir) config- 
uration is stable, while the (0, 7T, 7t) is at most mctastable. 
The DW can then be driven back and forth in the con- 
striction by an alternating current. That is an explicit 
realization of current-driven DW motion. 

We have then investigated the contribution of atomic 
relaxation to the DW migration barrier. The bottom 
panels of Figure show AW = W - W and AVFuni = 
W — Wuni as functions of 02 at different bias voltages V, 
where W(4>2) is the DW migration work (see Eq. 0) for 
atomic structure relaxed at the given bias and MS, while 
H^o(02) and Wu n i(02) correspond to atomic structure re- 
laxed at V = and to a uniform (unrclaxed) structure 
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J=leV J=1.5eV J=2eV 




FIG. 3: (Color online) DW-migration energy barriers at dif- 
ferent bias V for Jdd = 50 meV and J = 1, 1.5, 2 eV (panels 
from left to right). The solid (dashed) lines are for struc- 
ture relaxed at (0,0, n) [(0, 7r,7r)]. W(V,(f>2) corresponds to 
geometry relaxed at the given bias V; AWo = W — Wo and 
AWuni — W — Wuni, where Wo refers to geometry relaxed at 
V = and W un i to uniform geometry. 



J=leV J=2eV 




V(V) V(V) 

FIG. 4: (Color online) I — V curves for a geometry relaxed 
at V . Solid (dashed) lines represent (0, 0, 7r) [(0, 7T, 7r)] state. 
Top inset: dependence of net current for V — 1 V on the DW- 
migration reaction coordinate cj>%, circles represent structure, 
relaxed at <j>2 [see Fig. HJb)]. Bottom inset: A7 = I — io/uni, 
where Io (Iuni) refer to relaxed at V = (uniform) structure. 



respectively. The small current-induced atomic displace- 
ments (see Fig. 0) systematically increase the DWMB 
height (AWo > 0) by a few percent. AW does not show 
a dependence on the MS. 

However, the relaxation at V = from the uniform ar- 
rangement, which shortens the interatomic distance be- 
tween the live MMs by about 4% [Fig. (a)] , reduces 
the height of DWMB approximately by AW uni /W(4> 2 ) = 
25 - 30% (sec Fig. EJ for J = 1.5, 2 eV. This effect can 
increase dramatically (to about 200%) for exchange pa- 
rameters close to magnetic phase transitions. jjj, |l7J 
Thus the magnetic properties are strongly affected by 
the geometry of the contact, especially in the region of 
parameters where the J coupling mechanism starts com- 
peting with the direct exchange mechanism. 

Finally, the effect of the structural relaxation on the 
conductivity of the system is found to be small (Fig. 0}. 
The current is clearly insensitive to the DW migration 
within the constriction. The overall variation of the net 
current for the rotation of 02 for fixed geometry, which 
is in itself a small quantity, is further substantially com- 
pensated by structural rearrangement, i.e. the structure 
is found to respond to the magnetic rearrangement by 
structural adjustment, which minimizes the variation in 
the conductivity (see inset of Fig. Q). We have also 
found a decrease in conductivity due to relaxation of the 
structure from the uniform geom etry . This agrees quali- 
tatively with the findings in Rcf . |l5j , although the effect 
we observe is much smaller in magnitude. 



In conclusion we have developed a method to investi- 
gate the interplay between the magnetic and structural 
degrees of freedom of a ferromagnetic atomic point con- 
tact under bias. We have used it to assess the effects of 
the structural relaxation on the magnetic DWMB and 
reversely the effect of the magnetic configuration on the 
structural relaxation. Our main finding is that the in- 
terplay is only in one direction, that is the structural 
relaxation strongly modifies the DW migration barrier. 

In particular, we have found that the DWMB shows 
a substantial asymmetry, which increases with the ex- 
ternal bias even for a spatially symmetric system. That 
opens the possibility of voltage-controlled DW motion 
in such systems. The current-induced displacements 
(Ayy/a < 0.7%) from the relaxed at V — struc- 
ture produce a small positive shift in the DWMB height 
(~ 3%). This is small compared to the effect of the relax- 
ation from the initial uniform atomic configuration at the 
given bias. The latter is Ay un ;/a < 4% and has a much 
more dramatic effect on the DWMB profile, reducing the 
barrier height by upto 2/3 or even making the alterna- 
tive MS unstable, i.e. blocking the DW migration, for 
some exchange parameters. That is a signiture of the 
strong non-linear dependence of spin-polarized transport 
properties on structural rearrangements. However, struc- 
ture is not affected by the DW migration under bias. 
This is so because interatomic forces depend on the total 
charge density of the current-carrying electrons, not on 
their spin polarization. 
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